Tetra Point Wetting at the Free Surface of Liquid Ga-Bi 
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A continuous surface wetting transition, pinned to a solid/liquid/liquid/vapor tetra coexistence 
point, is studied by x-ray reflectivity in liquid Ga-Bi binary alloys. The short-range surface potential 
is determined from the measured temperature evolution of the wetting film. The thermal fluctuations 
are shown to be insufficient to induce a noticeable breakdown of mean-field behavior, expected in 
short-range-interacting systems due to their d u = 3 upper critical dimensionality. 
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The wetting transition, predicted independently by 
CahnQ and by Ebner and SaamQ in 1977, has at- 
tracted much theoretically and experimental]^] atten- 
tion, due to its importance for fundamental physics || 
and for applications || . Almost all experimental and 
theoretical studies of the wetting transition published to 
date address liquids interacting through long-range van 
der Waals potentials, for which the upper critical dimen- 
sionality is du < 3. Thus, for a d = 2 surface on a 
d = 3 bulk mean-field (MF) behavior is expected. By 
contrast, wetting phenomena in liquids interacting by a 
short range potential (SRW), where d u — 3, allow, in 
principle, to explore the regime where the MF behavior 
breaks down due to fluctuations and the renormalization- 
group (RG) approach becomes applicable. However, due 
to extremely demanding experimental requirements and 
the lack of a theory allowing to predict the SRW interac- 
tion parameters only a single experimental study address- 
ing fluctuation effects on SRW was published to date||. 
Even that study employs van-der-Waals-interacting liq- 
uids, and SRW was achieved by tuning the wetting tem- 
perature, TV, very close to the critical point, Tc- The 
authors' estimates of the fluctuations' magnitude pre- 
dicted a RG behavior. Nevertheless, a MF behavior was 
found and interpreted as an "apparent failure" of RG 
theory to account for SRW. 

We present here an x-ray study of wetting in a Ga-Bi 
binary alloy, which interacts by a short range screened 
Coulomb potential at all concentrations. This allows 
achieving SRW far from Tc, thus minimizing any pos- 
sible influence of the criticality on the wetting transition. 
Moreover, the x-ray techniques employed allowed the de- 
termination of the wetting layer's internal structure on 
an Angstrom scale and demonstrated the strong compo- 
sitional gradient within the wetting film, predicted by 
density functional calculations^, f7|. The wetting transi- 
tion in this system is pinned to a bulk monotectic tem- 
perature, Tm, creating a rare case of tetra point wetting, 
where four phase coexist in the bulk: a Ga-rich liquid, 
a Bi-rich liquid, a solid Bi and a vapor || |9|. We deter- 
mined the temperature variation of the wetting layer's 
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FIG. 1: The atomic fraction(c)-Temperature(T) bulk phase 
diagram of Ga-Bi. The symbols indicate coexistence lines of 
Predel's phase diagram lid]. The lines show the phase bound- 
aries calculated from thermodynamic data[23|. The dashed 
lines represent the metastable extension of the 1/1 coexistence 
line below Tm- The points are: C-bulk critical point, M- 
monotectic point, A,B,D-points on the experimental path. 
The insets illustrate the surface and bulk phases. In region 
II the wetting film is 50 A thick and the Ga-rich fluid is 5 
mm thick. The bold surface lines in the insets symbolize the 
Bi-monolayer. 



composition, thickness and internal structure along a tra- 
jectory in the phase diagram which includes several phase 
boundaries and the tetra point. The study highlights the 
intimate relation between the wetting layer's structure 
and the bulk phase diagram. The T-dependence of the 
wetting film's structure yields the first quantitative deter- 
mination of a MF short-range surface potential governing 
the wetting transition. Finally, from a RG analysis we 
estimate the influence of thermal fluctuations on the ob- 
served tetra point wetting. 

The measured bulk phase diagram of Ga-Bi is 
shown in Fig. |l|, with cartoons relating the surface phases 
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to the bulk ones. For all T a Gibbs-adsorbed Bi mono- 
layer is found at the free surface JO, Il2j. In Region I, 
T < Tm, the Bi monolayer and the solid Bi coexist with 
a Ga-rich liquid. In Region II, T M < T < T c , the 
bulk separates macroscopically into immiscible Ga-rich 
and Bi-rich liquids, and the denser Bi-rich phase sinks 
to the bottom of the sample pan. However, the Bi-rich 
phase also wets the free surface by intruding between the 
Ga-rich low density phase and the Bi monolayer, in de- 
fiance of gravity. In the homogeneous part of the phase 
diagram, region III, only the Bi monolayer is found at 
the free surface |L3|. The four phases, solid Bi, Bi-rich 
liquid, Ga-rich liquid, and vapor coexist at the boundary 
between regions I and II, Tm, rendering it a solid-liquid- 
liquid- vapor tetra point. 

X-ray reflectivity measurements were carried out at 
beamline X22B, NSLS, at a wavelength A = 1.54 A. The 
intensity R(q z ) reflected from the surface, is measured as 
a function of the normal component q z of the momentum 
transfer and allows determining the surface-normal elec- 
tron density profile p(;z)[fi4). The Ga-Bi alloy was pre- 
pared in an inert-gas box using > 99.9999% pure metals. 
A solid Bi pellet was covered by an amount of liquid Ga 
required for a nominal concentration c nom = 88 at% Ga. 
It was then transferred in air into an ultrahigh vacuum 
chamber. A day of bake-out yielded a pressure of 10 -10 
torr. The residual surface oxide on the liquid's surface 
was removed by sputtering with Ar + ions. Using ther- 
mocouple sensors and an active temperature control on 
the sample pan and its adjacent thermal shield a temper- 
ature stability and uniformity of ±0.05°C was achieved 
near Tm- 

X-ray reflectivity R(q z ) was measured at selected tem- 



peratures on the path A— >M. The standard procedure) 15 



for determining p(z) from the measured R(q z ) is to con- 
struct a physically motivated model for p(z) and fit 
its Fourier transform to the measured R(q z ) / 'Rf{<1z)- 
Rp{q z ) is the Fresnel reflectivity from an ideally flat and 
sharp surface having the electron density of the Ga-rich 
liquid. We employ a two-box model) 11 , where the upper 
box represents the Gibbs adsorbed Bi monolayer, and 
the lower box - the Bi-rich wetting film. The model 
also includes three roughnesses for the three interfaces: 
vapor/Bi- monolayer, Bi-monolayer/Bi-rich film, and Bi- 
rich film/Ga-rich subphase. The fits (lines) to the mea- 
sured (points) R/Rf are shown in Fig. ||(a), and the 
corresponding p(z) profiles - in Fig. ||(b). At point A 
(T A = 205°C, inset to Fig.|(a)), typical of region I, R/R F 
exhibits a pronounced maximum at q z — 0.8A -1 , in- 
dicative of a high electron density layer at the surface. 
The p(z) obtained from the fit is consistent with a seg- 
regated monolayer of pure Bi at the free surface [ll] fl2| . 
As the temperature is increased towards M, a peak in 
R/Rf starts to develop at q z « 0.2A -1 , gradually shift- 
ing to lower q z . This behavior manifests the formation 
of the wetting layer and its continuous growth in thick- 
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FIG. 2: (a) R/R F for selected temperatures T A < T < Tm- 
(♦) 205.0°C = T A , (<) 218.9°C, (□) 220.4°C, (0) 221.0° C, 
(V) 221. 5°C, (A) 221.8° C, (•) 222.0°C = T M - Solid lines: 
fits to a two-box model of p(z). (b) Fit-refined electron den- 
sity profiles p(z) for aforementioned T's with Ta < T < Tm, 
Psub = electron density of Ga-rich subphase. A and M mark 
p(Ta) and p(Tm), respectively. 



ness upon approaching Tm- For the fully- formed wetting 
layer at M (Tm = 222.0°C), R/R F shows two \ow-q z 
peaks, characteristic of a thick film, and the fit yields a 
maximal film density of p I p sub-phase. = 1-26 as shown in 
Fig. ||(b). This agrees well with the value of 1.25 calcu- 
lated for this ratio from the phase diagram at Tm- The 
fit-refined ~ 50A film thickness also agrees with ellip- 
sometric measurements fl6f . As can be seen in Fig. ||(b), 
the evolution of the density profiles along the path A^M 
proceeds via highly inhomogcneous wetting films that are 
similar to the microscopic density profiles calculated us- 
ing density functional theory [|[ 

The surface and bulk transitions at Tm are driven by 
the excess free energy, A/i m , of the metastable Bi-rich 
phase over that of the Ga-rich phase. Thus, to relate 
the observed continuous surface transition to the first- 
order bulk transition at Tm, we replot in Fig. |](a) the 
(c,T) phase diagram of Fig. [I] on the (Ap m ,T) plane. 
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FIG. 3: (a) (A/i m ,T) phase diagram: (A-M) is the liquid- 
solid, and (M-C) is the liquid/liquid coexistence lines. The 
path B-D is in the single phase region, and M and C are the 
monotectic and critical points. Inset: effective wetting layer 
thickness d on A^M (squares) and B— >D (open circles). The 
solid line is a fit to the A— >M d values, (b) The measured d 
along the experimental path. 



For T > Tftf, the 1/1 coexistence line transforms into a 
straight line at (Ap m = 0,T) extending from M to C. 
For T < Tm the metastable 1/1 coexistence line (dashed 
line) extends horizontally below Tm, while the solid- 
Bi/Ga-rich coexistence line goes below the metastable 
line from M to A. This illustrates Dietrich and Schick's 
observation[|| that the path leading to coexistence on 
heating from A to M is dictated by the topology of the 
phase diagram and corresponds to a path which probes 
"complete wetting": Ap m — > on the path A — > M. 

To determine the surface potential, we define an 
effective wetting film thickness, d = f°°[p(z) — 
PGa-nch\/[pBi-rich - PGa-rtch]dz. Here z s is the top of 
the wetting film as marked in Fig. ^(b), and psi-rich, 
PGa-rich are the electron densities of the coexisting bulk 
liquid phases, calculated from the phase diagram. The d 
values calculated from the reflectivity fits along A— >M, 
M->B, and B^D are plotted in Fig. |(b). Along A^M 
a continuous increase in d is observed, while on the on- 
coexistence M— >B path a constant d w 50 A is found. 
This is in agreement with the predictions for SRW: a con- 
tinuous d growth limited by gravity effects B. At B the 
Ga concentration reaches its nominal value, c nom = 88 



at%, in the homogeneous liquid phase (region III). Upon 
further increase of T c remains unchanged. The sample 
leaves the 1/1 coexistence line, moving further into region 
III, and the wetting film vanishes continuously. Further 
details of the M— >D path are discussed elsewhere p|. 

As implied by the similarity to Ebner and Saam'sQ] 
results, a detailed calculation of the wetting film's inho- 
mogeneity will require either a density functional anal- 
ysis, or some other equivalent approach. Nevertheless, 
even a simple model approximating the wetting layer by 
a slab of thickness d allows a confident determination of 
the surface potential. For this model, the grand canonical 
potential tts at the surface is given in the MF approxima- 
tion by n s = NA [dAp + £,MF®e~ d/ Z MF ] with a short- 
range potential of decay length £mf and amplitude <E>. 
Here Ap = Ap m +Ap g is the excess Gibbs Free Energy of 
the Bi-rich wetting phase over that of the Ga-rich phase, 
and Ap g — g Ap h V m is the gravitational energy paid 
for having the heavier Bi-rich phase at the surface. The 
particle number density per unit volume is N, Aq is an 
arbitrary surface area, Ap is the mass density difference 
of the two phases, V m is the molar volume of the Bi-rich 
phase, and h is the height difference between top and bot- 
tom. The equilibrium thickness, d* , is that minimizing 
fig'- ^fU* = 11 ■ Neglecting the gravitational term for 
points off 1/1 coexistence yields: d* = £mf m ($/A^i m ). 
The inset in Fig. ||(a) shows the experimental d values on 
the path A^M (solid squares). The clear linear depen- 
dence on ln(l/A/i) confirms the major theoretical pre- 
dictions for SRW: a continuous growth and a logarith- 
mic divergence of the wetting film's thickness A fit 
of d* (line) to the experimental d (solid points) yields 
£,mf = 6.3 A and $ = 43 J/mol. Both values are rea- 
sonable for such a metallic system. The value of <E> cor- 
responds to a change in surface energy upon wetting of 
about 400 mJrn -2 which is consistent with surface en- 
ergy measurements on Ga-Bi alloys [|l7), and Ga-Pb ones 
H . It is interesting to compare the behavior of d for the 
path A— >M with that along the path B^D, shown in 
the inset of Fig. || (open circles). The values overlap for 
most of the path, demonstrating that the behavior along 
A— >M is the same as can be observed for any other path 
probing complete wetting, e.g. the path B— >D which is 
solely determined by the choice of c nom . 

The analysis above exhibits a very good agreement of 
the experimental SRW results with the universal behav- 
ior predicted by MF theory. However, it is important to 
note that this does not automatically mean a failure of 
the RG predictions for a non-universal SRW behavior in 
our, and other, d u — 3 systems. Both MF and RG ap- 
proaches yield here a continuous, logarithmic divergence 
of d* at Tyy, but with different amplitudes for p"8"| . 
For MF, the amplitude is £mf while for RG, where d* ~ 
£ RG (I + oj/2) \n(l/Ap), the amplitude is £ RG (l+w/2). 
The "capillary parameter" u> = fcgIV/(47r£^7) (where 
£b is the bulk correlation length, and 7 is the 1/1 in- 
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terfacial tension) measures the magnitude of the dom- 
inant thermal fluctuations in the system, the thermally 
induced capillary waves at the 1/1 interface of the coex- 
isting Bi- and Ga-rich liquids. Hence, w measures the 
deviation from a MF behavior, which, of course, neglects 
all fluctuations. A calculated 7 = AmJmT 2 and a 
£b = 6.1 ± 1 A estimated from the two-scale factor uni- 
versality (TSFU) theory for bulk demixing in Ga-Bijl9[, 
yield w(T M ) = 0.3 ± 0.2. This is well within the cu < 2 
range, where RG analysis predicts only a mild effect on 
complete wetting Q. Indeed, £rg = 5.4 A is only ~ 15 
% smaller than £mf — 6.3 A above. The two values are 
close, and in reasonable agreement with the value of 
expected from the TSFU||. Thus, a clear distinction 
between RG and MF behavior can not be drawn in our 
case. 

Finally, we compare the measured gravity-limited max- 
imal thickness of 50 A along M— >B with theory. On coex- 
istence A/i m = and thus d* = £rg hi ($/A/x s ). Using 
$ = 43 J/mol and the known material constants in A/i g 
yields d* rav = 15.6 ^rg ~ 85 A. Since this analysis does 
not take into account the excess energy associated with 
concentration gradients across the interfaces some over- 
estimation of d* is not too surprising. However, even this 
rough calculation shows that the wetting film thickness 
is expected to be on a mesoscopic rather than on the 
macroscopic length scale observed for similar wetting ge- 
ometries in systems governed by long-range, dispersion 
forces [[|. 

We presented here the first detailed, Angstrom resolu- 
tion study of a short-range wetting transition at the free 
surface of the binary Ga-Bi liquid alloy. The character- 
istic properties of SRW, a logarithmic divergence of the 
wetting layer's thickness upon approaching the transi- 
tion, are clearly exhibited by the measurements. The ob- 
served strongly inhomogeneous internal structure of the 
wetting layer is similar to that obtained from density 
functional calculations of wetting at a hard wall|^|, 0. 
The surface wetting transition is found to be pinned to 
the bulk monotectic point M in this alloy, and the topol- 
ogy of the bulk phase diagram enforces a path probing 
complete wetting^]. Since this path ends in M where 
four phase coexist, the surface transition constitutes a 
rare tetra point wetting effect. An estimation of the 
capillary parameter u> indicates that fluctuations affect 
the observed complete SRW transition only marginally. 
A critical wetting transition with a similar value of u> 
should show more pronounced deviations from MF be- 
havior. Here, however, this critical transition occurs on 
the metastable extension of the 1/1 coexistence line some- 
where below Tm, and may be reachable only if the al- 
loy supercools sufficiently. Due to the mild effect of the 
fluctuations, the surface potential driving the transition 
could be determined quantitatively. This is of particular 
importance since in contrast with van-der-Waals interact- 



ing liquids, where the surface potential can be estimated 
from the Dzyaloshinski-Lifshitz-Pitaevskii theory^, pl| , 
no equivalent theory is available for liquids with short- 
range interaction. We hope that this first quantitative 
determination of the SRW potential will stimulate the 
development of such theories. Available predictions, in 
turn, would initiate more experimental studies of fluc- 
tuation effects in SRW. In particular, they may allow 
a quantitative investigation of Parry's prediction of an 
additional renormalization of lu in SRW due to a cou- 
pling between the fluctuations at the liquid/vapor and 
at the 1/1 interfaces ]2(J. This coupling may resolve the 
dichotomy between RG prediction and the observations 
in the single previous experimental study of SRW|(|. 
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